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Abstract 

Background: Prostate cancer is characterized by heterogeneity in the clinical course that often 
does not correlate with morphologic features of the tumor. Metastasis reflects the most adverse 
outcome of prostate cancer, and to date there are no reliable morphologic features or serum 
biomarkers that can reliably predict which patients are at higher risk of developing metastatic 
disease. Understanding the differences in the biology of metastatic and organ confined primary 
tumors is essential for developing new prognostic markers and therapeutic targets. 
Methods: Using Affymetrix oligonucleotide arrays, we analyzed gene expression profiles of 24 
androgen-ablation resistant metastatic samples obtained from 4 patients and a previously published 
dataset of 64 primary prostate tumor samples. Differential gene expression was analyzed after 
removing potentially uninformative stromal genes, addressing the differences in cellular content 
between primary and metastatic tumors. 

Results: The metastatic samples are highly heterogenous in expression; however, differential 
expression analysis shows that 4 1 5 genes are upregulated and 364 genes are downregulated at least 
2 fold in every patient with metastasis. The expression profile of metastatic samples reveals changes 
in expression of a unique set of genes representing both the androgen ablation related pathways 
and other metastasis related gene networks such as cell adhesion, bone remodelling and cell cycle. 
The differentially expressed genes include metabolic enzymes, transcription factors such as 
Forkhead Box Ml (FoxMI) and cell adhesion molecules such as Osteopontin (SPPI). 
Conclusion: We hypothesize that these genes have a role in the biology of metastatic disease and 
that they represent potential therapeutic targets for prostate cancer. 
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Background 

Prostate cancer is the most common cancer in men result- 
ing in over 232,090 new cases and 30,350 deaths annually 
[ 1 ] . For prostate cancer patients, metastatic disease reflects 
the most adverse clinical outcome. Osseous involvement 
with severe bone pain and spinal cord complications 
occur commonly in patients with metastatic disease [2J. 
However there is considerable heterogeneity in outcome 
after primary diagnosis and currently there are no mor- 
phologic or circulating biomarkers that can accurately pre- 
dict the development of metastatic disease. 

Metastatic prostate cancer represents the tumor's ability to 
escape from the primary organ and eventually colonize a 
distant site. Disruption of a complex set of biological 
processes must occur in order for tumor cells to leave the 
prostate and establish themselves in a different environ- 
ment. Their altered interaction with the prostate microen- 
vironment, including the stroma and extracellular matrix, 
their ability to migrate into the vasculature and establish 
themselves in secondary organs with recruitment of vascu- 
lar supply represent disruption of normal cellular proc- 
esses [3]. Understanding the molecular events involved in 
the development of metastatic prostate cancer has the 
potential to identify biological determinants that can aid 
in prognosis and development of more effective therapies. 

Using gene expression microarrays, a number of studies 
riave characterized expression profiles of prostate cancer, 
riormal tissue and metastatic cancers. In some cases, cor- 
relations between tumor expression signatures, clinical 
parameters and outcome have been identified [4-11]. 
Unique profiles have been reported for untreated and 
short-term androgen ablation treated organ-confined dis- 
ease and for metastatic disease, with a subset of genes dif- 
ferentiating metastatic androgen ablation resistant 
prostate cancer (AARPC) from androgen dependent met- 
astatic cancers [10,12-14]. In general, metastatic prostate 
cancer is characterized by changes in expression of genes 
involved in signal ttansduction, cell cycle, cell adhesion, 
migration and mitosis. In addition to these genes, AAR- 
PCs exhibit changes in expression of the androgen recep- 
tor and enzymes involved in the sterol biosynthesis 
pathway [12]. 

Some of the genes previously reported as highly downreg- 
ulated in prostate tumors may reflect the differences in 
cellular content of metastatic and organ-confined tissues 
rather than intrinsic differences in biology. In contrast 
with organ-confined prostate tumors which are composed 
of a mixture of glandular epithelial, smooth muscle and 
other stromal cells, metastatic tissue samples are almost 
exclusively epithelial, with minimal supporting stroma 
and absence of smooth muscle. In this study, we charac- 
terize gene expression in androgen ablation resistant met- 



astatic tumors after removing potentially uninformative 
suomal genes. The deleted sttomal genes consist of those 
reported in a recent report characterizing the gene expres- 
sion patterns in the prostate stroma, tumor and normal 
epithelium [15]. Our results provide novel insights into 
the biology of metastasis. 

Methods 

Tumor sample procurement 

All tissue samples were acquired from the Health Sciences 
Tissue Bank of the University of Pittsburgh Medical Center 
under stringent Institutional Review Board guidelines 
with appropriate informed consent. The 18 donor and 64 
primary prostate tumor samples have been described pre- 
viously [7]. Specimens were received directly from the 
operating room. Samples (>500 mg) were excised and 
snap frozen in liquid nitrogen within 30 min of excision 
and stored at -80 °C until extraction of RNA. Metastatic 
tumor samples were obtained from a warm autopsy pro- 
gram and processed similarly to primary tumors. An H&E 
stained frozen section of each sample was evaluated by a 
pathologist, to determine epithelial and stromal content 
and verify the presence of tumor in the sample. Dissection 
of the frozen tissue block was performed with the guid- 
ance of a marked H&E slide to minimize the presence of 
host tissue in the metastatic samples. All samples used in 
the study contained >80% tumor. Metastatic tumor sam- 
ples were minced and divided into two equal portions to 
be extracted with the sample protocol used for each set of 
primary tumors. 

Clinical profile of cases 

The clinical characteristics of the 64 primary tumor sam- 
ples used in the Affymetrix portion of our study have been 
previously described [6,7]. These cases have a mean fol- 
low-up time of 3 years. The metastatic samples consisted 
of 24 tissues derived from 4 patients (Table 1). All patients 
with metastatic disease had received androgen ablation 
therapy and had shown progression of disease while on 
androgen ablation. The clinical characteristics for the 
additional 10 primary prostate tumor cases used in the 
CodeLink study are shown in Table 1. 

RNA extraction 

RNA purification for the 64 primary samples has been pre- 
viously described [6]. The set of metastatic samples ana- 
lyzed with the Affymetrix platform was extracted with the 
same methodology. The set of metastatic samples and pri- 
mary tumors analyzed with the CodeLink platform were 
exteacted using the RNeasy kit (Qiagen, San Diego, CA). 
For the metastatic samples, one sample did not have 
enough for extraction with the Qiagen method, only 23 
metastatic samples are included in the CodeLink assays. 
The concenttation of each total RNA sample was meas- 
ured with a Nanodrop ND-1000 spectrophotometer 
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fable I: Clinical variables for primary and metastatic prostate ce 



;r samples used in this study 



f Prostate Cancer 
'■' Tissue Samples 


No. of Samples 


Number of 
Patients 


Microarray 
Platform 


Clinical Information 




64 


64 


Affymetrix 


Please see re 1 


ference [7] 










Gleason Score 


No. of Cases 










7 


10 


Primary Tumors 


10 


10 


CodeLink 


Pathological Stage 


No. of Cases 










2B 
3A 
3B 


6 
3 










Patient ID 


No. of Samples 










FB656I 
FB666 
FB667 
FB669 


1 

8 
5 



Para Aortic Lymph 

Node 
Para Tracheal Lymph 
Node 



Lymph Node 



(Nanodrop Technologies, Wilmington, DE). RNA integ- 
rity was determined by capillary electrophoresis using an 
Agilent 2100 Bioanalyzer (Agilent, Willmington, DE). 

cRNA preparation and gene expression assays 

cRNA was prepared and hybridized to Affymetrix Gene- 
Chip HGU95av2, HGU95b and HGU95c arrays (Affyme- 
trix, Santa Clara, CA) as previously described [6]. For gene 
expression profiling with the CodeLink Gene Expression 
System (GE Healthcare, Piscataway, NT), biotin-labeled 
cRNA was prepared as previously described [16]. Ten 
micrograms of biotin-labeled cRNA product from each 
sample were then fragmented with RNA fragmentation 
buffer at 94 0 C for 20 minutes. Hybridization mix was pre- 
pared according to the manufacturer's instructions and 
the final volume was adjusted to 260ul using nuclease- 
free water. The hybridization mix was heat denatured at 
90 °C for 5 minutes, cooled on ice and then applied to 
Human Uniset 20 K arrays (GE Healthcare, Piscataway, 
NJ). Arrays were incubated at 37° C for 18 h with shaking 
at 60 rpm in an Innova hybridization oven (New Bruns- 
wick, Edison, PA). 



After hybridization, arrays were placed in a pre-heated 
(46 • C) chamber filled with 0.75 x TNT (0.75 M Tris-HCL, 
pH 7.6, 3.75 M NaCl, Tween-20, and milli-Q water) and 
incubated at 46° C for 1 hour. Arrays were then stained 
with Streptavidin-Alexa Fluor 647 (Molecular Probes, 
Grand Island, NY) for 30 minutes at room temperature. 
Upon the completion of staining, the arrays were washed 
three sequential rimes in fresh 1 x TNT (1 M Tris-HCL, pH 
7.6, 5 M NaCl, Tween-20, and milli-Q water) and then 
washed two final times in fresh solutions of 0.05% 
Tween-20 and 0.1 x SSC with gende agitation. All arrays 
were dried by centrifugation at 2000 rpm for 3 minutes. 

Affymetrix arrays were scanned in an Affymetrix GCS3000 
Scanner (Affymetrix, Santa Clara, CA). CodeLink arrays 
were scanned with the GenePix 4000B scanner using 
GenePix Pro 4.1 software (Molecular Devices, Sunnyvale, 
CA). 

Gene expression data analysis 

The raw scanned array images from the Affymetrix Gene- 
Chip U95 arrays were processed using GCOS 1.1 software 
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rising the MAS5 algorithm (Affymetrix Corporation, Santa 
Clara, Ca) to generate probe eel intensity (*.cel) files. Data 
normalization to remove variation in overall chip intensi- 
ties was performed by global scaling to a chip mean target 
intensity of 200 (MAS 5.0). Data for U95Av2, B and C 
arrays were combined for further analyses. 

To identify differentially regulated genes in both datasets, 
these were analyzed with the Significance Analysis of 
lilicroarrays software (SAM v 1.2) [17]. Prior to analysis, 
genes that showed low variation across all samples were 
removed by using the filtering option in the Avadis 3.3 
Pride Software (Strand Life Sciences, Bangalore, India) 
data analysis tool. To avoid false results due to difference 
in the tissue composition of metastatic and primary 
tumors, genes identified as being highly expressed in the 
prostatic stroma as per Stuart etai [15] were also removed. 
In all 1506 stromal genes and 7678 invariant genes were 
removed from the Affymetrix dataset. SAM generated gene 
lists with the lowest false discovery rates (FDR) were fur- 
ther analyzed for gene ontology (GO) and pathway anno- 
tations using NIH's DAVID annotation tool [18]. 

For CodeLink arrays, image files were analyzed with the 
GodeLink Expression Analysis Software version 4.1 (GE 
Healthcare) with use of the normalized intensity values 
in downstream analysis. For cross-platform comparison, 
Affymetrix probe sets and Codelink identifiers were 
inapped to Unigene ids using the DAVID annotation tool 
(see above). Expression data from both platforms was 
compared using z-transformation. Hierarchical clustering 
vfas performed using Eisen's Cluster and Treeview [19]. 
Data from Affymetrix experiments has been submitted to 
NCBI's Gene Expression Omnibus (GEO) as series 
GSE6919, with the following accession numbers 
GSE6604 (normal donor prostate), GSE6605 (metastatic 
prostate tumors), GSE6606 (primary prostate tumors) 
and GSE6608 (normal prostate tissue adjacent to tumor). 
Data from the CodeLink platform have been submitted to 
GEO with the accession number GSE6752 (primary and 
metastatic prostate tumors). 

Quantitative real-time PCR 

Differential expression of ten genes in primary and meta- 
static prostate cancer samples was verified with quantita- 
tive real-time PCR (QPCR) with the ABI PRISM 8 7000 
sequence detection system (Applied Biosystems, Foster 
City, CA). Three selected RNA samples from each patient 
Were pooled together (except for patient FB666 n = 1) and 
therefore four RNA samples, each representing one 
patient, were tested. RNA samples were first heat-dena- 
tured at 70 °C for 10 minutes in an Eppendorf master 
cycler (Eppendorf, Westbury, NY) and then chilled imme- 
diately on ice. cDNAs were reversely transcribed from one 
microgram of RNA using the M-MLV Reverse Tran- 



scriptase kit (Invitrogen, Carlsbad, CA), as recommended 
by the manufacturer. QPCR was performed based on the 
manufacturer's instructions with TaqMan Gene Expres- 
sion Assays (Applied Biosystems) for the following genes: 
EGR3, SYNP02, ANGPT2, SPP1, FOXM1, ADM, RDX, 
TGFBRAP1, MAK and EGRl(assay IDs: Hs00231780_ml, 
Hs00326493_ml, Hs01048047_mH, Hs00959010_ml, 
Hs00153543_ml, Hs00969450_gl, Hs00988414_gl; 
Hs01093285_ml; Hs01048300_ml, Hs00152928_ml). 
When multiple TaqMan assays for one gene were availa- 
ble, the assay that interrogated the sequence closest to the 
target sequence in the Affymetrix arrays was chosen. PCR 
cycles were performed according to the assay instructions 
in an ABI PRISM* 7000 Sequence Detection System 
(Applied Biosystems). Relative quantification of the 
expression level of each transcript in each sample was cal- 
culated using the Delta-Delta CT method in the ABI 
PRISM 7000 Sequence Detection System Software 
(Applied Biosystems) [20]. Human reference RNA from 
Suatagene (Stratagene Corp., La Jolla, CA) was used as the 
calibrator (untreated control) and human glucuronidase 
beta (GLTSB) gene was used as the endogenous reference 
gene (Forward primer: GGAAITITG CCG ATTTCATGA; 
Reverse primer: CCG AGT GAA GAT CCC CTTTTT; Probe: 
6FAM-AAC AGT CAC CGA CGA GAG TGC TGG G- 
TAMRA). 

Results 

Differential gene expression in metastatic prostate cancer 
and the role of stromal content in defining true 
downregulated genes 

Differential expression analysis of the metastatic and pri- 
mary tumor samples shows that a large number of the 
most highly downregulated genes such as TAGLN, 
ACTG2, TPM1, MYH111 and DES have been previously 
identified as expressed mostly in the prostatic stromal 
cells [15]. Since only the epithelial component of prostate 
cancer is present in metastatic tumors, this result most 
likely reflects the lack of stroma in metastases, and not a 
true down-regulation of these genes in the metastatic epi- 
thelial cells. Therefore, based on a recent report character- 
izing cell type specific gene expression in the prostate [ 15], 
we removed the set of genes expressed mainly by the stro- 
mal cells of the primary tumors. In all 1506 transcripts 
associated with a stromal signature were deleted prior to 
further analysis. Since the stromal genes were character- 
ized using the U95Av2 chip and our analysis includes 
u95Av2, B and C chips, only stromal genes represented by 
probe sets on U95Av2 were removed in this modified 
analysis. SAM analysis shows that 1277 genes are up and 
977 genes are downregulated at least 2 fold at the lowest 
FDR (0.01), in metastatic prostate samples (see Addi- 
tional file 1). A list of the top 50 up and top 50 down reg- 
ulated genes at the lowest fdr, after removing ESTs and 
uncharacterized clones is shown in Table 2. This list 
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includes signal transducers, cell cycle regulators, meta- 
bolic enzymes and cell adhesion molecules. Some of the 
most upregulated genes in our list are EIF1AX, AR, HSPD1 
arid HSPCA, K-ALPHA1, MLL5, UGT2B15, and some of 
the most downregulated genes include WNT5B5, 
ANXA1 1, POS and SPRP1 . 

Metastatic samples are heterogenous in gene expression 

Using immunohistochemistry (IHC) Shah et al. have 
shown that metastatic samples are highly heterogenous in 
expression of prostate specific markers leading to the 
hypothesis that at the molecular level, metastatic prostate 
cancer may represent multiple diseases even within the 
same patient [21]. We examined the expression of several 
transcripts markers including some studied by Shah et al. 
and confirmed the heterogeneity of expression levels in 
metastatic prostate cancer tissues. Expression values in 
donor samples, primary and metastatic samples were 
compared. Prostate specific antigen (PSA/KLK3) remains 
high in some metastatic samples and is low or absent in 
others, even within the same patient (Figure 1). Interest- 
ingly, AMACR, another biomarker for prostate cancer [22] 
expresses a heterogenous expression pattern similar to 
PSA. HPN, which is overexpressed in primary cancer 
maintains high expression in the metastatic samples in 
our study. AR, while overexpressed in 23, out of the 24 
metastatic samples, shows highly variable expression val- 
ues in individual samples. The proto-oncogenes FOS and 
j'tJNB, which are both overexpressed in primary tumors, 
are consistently downregulated in all metastatic samples. 

Genes regulated in all metastatic cases 

Hierarchical clustering analysis reveals that gene expres- 
sion in metastatic samples is more variable between 
patients than between different metastatic, sites from each 
patient (Figure 2). Although the 24 metastatic samples 
represent tissues from 6 metastatic sites (Table 1), no 
drgan specific clusters were detected (Figure 2) whereas 
samples from the same patient tend to cluster together. 
Statistical comparison of organ-specific expression pro- 
files was not attempted due to unequal distribution of 
samples from different metastatic sites. 

In order to identify probe sets that are similarly regulated 
in every patient, and therefore likely to represent a specific 
metastatic profile, the SAM differentially expressed gene 
list at a FDR of 2% was further filtered. For each gene on 
tins list, a patient specific median expression value was 
calculated from the multiple samples from each patient. 
Patient P4 had only sample and this sample's signal value 
was considered the median value. The median values were 
then compared to the median value of the primary sam- 
ples and those probe sets whose median value showed 
equal or more than a 2 fold change in every patient were 
considered part of the metastatic prostate cancer signa- 



ture. Under this criteria 415 transcripts are upregulated 
and 364 are downregulated in all patients with metastasis 
(see Additional file 2). A truncated gene list consisting of 
genes regulated at least 3 fold in all patients is shown in 
Table 3. Upregulation of AR in all samples from meta- 
static cancer patients represents a known " androgen resist- 
ant" or AARPC (androgen ablation resistant prostate 
cancer) phenotype [12]. The transcripts identified as dif- 
ferentially expressed in our study exhibit similarities with 
a previous study of AARPC tumors [10,12]. Cytokeratins 
5 and 15 [KRT5/KRT15), markers of basal cells in prostate 
glands, show uniform downregulation in all metastatic 
tumors, confirming the absence of basal epithelial cells. 

Biological annotation of differentially expressed genes in 
metastatic prostate cancer 

The list of differentially expressed transcripts at least 2 fold 
in all patients was further analyzed for biological themes 
and gene ontology (GO) using the NIH's DAVID annota- 
tion tool. This analysis revealed that metastatic prostate 
cancer exhibits altered regulation of amino acid, carbohy- 
drate and nucleotide metabolism consistent with the pro- 
liferative capacity and altered energy needs of metastatic 
tumors (data not shown) . In the context of prostate cancer 
biology, genes involved in cell-adhesion, bone remode- 
ling, cell-cycle and transcription are of particular interest 
(Table 4). Disruption of cell adhesion and altered interac- 
tion with the extracellular matrix is a hallmark of meta- 
static tumors [3]. In agreement with this, the secreted 
phosphoprotein and cell adhesion molecule osteopontin 
(SPP1) is one of the most highly upregulated transcripts 
in our metastatic samples. Elevated expression of SPP1 
has been correlated with poor prognosis in prostate 
tumors and other cancers and it has often been implicated 
in metastasis to bone and other organs [2,23-28]. In all 29 
probe sets representing cell adhesion genes are altered in 
all metastatic samples. This gene list includes FN1, ITGB8, 
THBS2, HNT and CDH10. Genes involved in bone 
remodeling such as BMP4 and ANKH are also altered in 
expression, although none of the samples in our study are 
bone metastatic samples suggesting that these proteins 
may also be involved in cancer metastasis to other organs. 

Disruption of the cell cycle is highlighted by the presence 
of a large number of cell-cycle related transcripts in the list 
of differentially expresed genes in all metastatic samples. 
The list contains 37 cell cycle genes, and includes SEP4, 
SEP7, PTN and VEGF. Similarly, a large number of tran- 
scription factors (67) including AR, SRY, FOS and EGR3, 
are differentially expressed. Two members of the winged- 
helix family of transcription factors, FoxPl and FoxtAl, 
show upregulation in the metastatic samples. Interest- 
ingly, FoxMlb has been shown to promote progression of 
prostate carcinomas in an experimental model [29]. 
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fable 2: Top 100 genes differentially expressed in metastatic samples compared to primary tumor samples 
Gene Symbol ProbeJD d_Value Fold Change 
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S60S6_at 


11.31 


2.64 
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SDCCAG3 


430l4_at 


9.91 


6.04 
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table 2: Top 100 genes differentially expressed in metastatic samples compared to primary tumor samples (Continued) 
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-3.23 


MRC2 
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-7.13 


-2.38 
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FLJ20542 


50l64_at 


-7.36 


-3.13 
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56272_at 


-7.42 
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-3.70 


FOS 
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FOS 
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PILB 
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59669_at 


-7.77 


-2.70 


RPS27L 


564IO_at 


-7.80 


-2.27 


t"OMM20 


36l98_at 


-7.91 


-2.13 


L 


438l9_£_at 


-7.91 


-3.13 


5TAT6 


4l222_at 


-7.92 


-2.78 


SELM 


64449_at 


-7.95 


-9.09 


BOC 


52999_at 


-8.15 


-4.17 


L 


44746_at 


-8.20 


-2.13 


SMBP 


46307_at 


-8.28 


-3.23 


'iiVNTSB 


6l330_at 


-8.34 


-3.33 


FLJ22386 


50l98_at 


-8.43 


-4.55 


^VNT5B 


58787_at 


-8.61 


-2.27 


^DHHC4 


45807_at 


-8.81 


-2.13 










kpiP 


38063_at 


-9.04 


-2.33 


ANXAI 1 


55664_at 


-9.16 


-2.94 


fvNT5B 


66l42_s_at 


-9.33 


-4.17 


^VAS 


38963_i_at 


-9.43 


-2.50 


JFCI 


44820_f_at 


-9.45 


-3.33 


JFC1 


48805_f_at 


-9.48 


-2.44 


WNTSB 


6l292_s_at 


-9.86 


-7.14 


GIRBP 


39864_at 


-10.94 


-3.70 



Gene expression data from the Affymetrix platform for 25 metastatic and 64 primary tumor samples was analyzed for differential gene expression 
fay SAM. The differentially expressed genes with the lowest FDR were sorted by fold change. The top 1 00 genes, organized by fold change are 
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Figure 2 

Hierarchical clustering of primary and metastatic prostate cancer samples. The 24 metastatic (Mets PI, Mets P2, 
Mets P3 and Mets P4) and 64 primary tumor samples were clustered. The top row of color coded boxes represents metastatic 
or primary samples; the bottom row represents the organ from which the sample was obtained. 



The MAP kinase signaling pathway was also identified as 
being important in the metastatic process; with 26 probe 
sets involved in this pathway being differentially 
expressed in all metastatic samples. The regulated genes 
include DUSP1, DUSP2, DUSP8, MAP3K8, MAP4K4, 
FGF13, FGFR2 and FOS. Involvement of MAP kinase in 
androgen receptor signaling has been previously 
described [30]. 

Validation of differentially expressed transcripts with an 
independent set of primary tumors and different gene 
expression platforms confirms gene expression profiles of 
metastatic prostate cancer 

Gene expression analysis with the CodeLink Uniset 20 K 
riiicroarray was carried out for 23 of the metastatic sam- 
ples and compared to an independent set of 10 primary 
ttrmors. Similar to the Affymetrix analysis (see above), 
hierarchical clustering of the CodeLink data set reveals 
heterogeneity in expression and no organ-specific cluster- 
ing (data not shown). Comparison of results with the 
Affymetrix based dataset, based on genes with common 
Unigene ids on both platforms, show a similar pattern of 
differentially expressed genes. Of the top 1000 up and 
down regulated transcripts from each platform, approxi- 
mately 70% share common unigene ids and of these 22% 
of the genes are identified as regulated by both platforms 
($iee Additional file 3). This level of correlation is signifi- 
cant, given the well-documented difficulties in cross-plat- 
form comparisons of expression data [31,32]. Examples 
of z-transformed expression values for selected genes in 
both platforms are shown in Figure 3. 

Additionally, real-time quantitative PCR .(QPCR) assays 
vfere performed for a selected set of genes in pooled sam- 
ples for each patient with metastatic disease and 5 of the 
primary tumors from the CodeLink set. The transcripts for 
this analysis were chosen to represent diverse biological 
processes and were chosen from the differentially 



expressed genes identified as up/down-regulated in the 
Affymetrix/Codelink data comparison. As shown in Fig- 
ure 4, qPCR assays confirmed the results from the micro- 
array platforms. SYNP02 and EGR3, which are 
downregulated and RDX and FOXM1, which are upregu- 
lated in the microarray analysis exhibit a very similar 
expression pattern in the qPCR analysis. Interestingly, 
FoxMl is consistently up regulated in metastases, while 
RDX was upregulated in only two of the four patients with 
metastatic disease, confirming the heterogeneity of meta- 
static prostate cancer. 

Discussion 

Despite extensive research, the molecular mechanisms of 
metastatic prostate cancer and androgen resistance devel- 
opment are still poorly understood. Our study shows that 
a number of biological processes including cell adhesion, 
cell cycle and transcription regulation are altered in meta- 
static disease when compared to primary tumors, and 
point to specific transcripts that participate in the meta- 
static process. 

Previous investigators have reported differences in gene 
expression profiles of metastatic and primary prostate 
cancer [10,12,14,21]. Our results show partial overlap 
with these previous characterizations of metastatic dis- 
ease. Some genes that are in concordance with these stud- 
ies include transcription factors such as FOXM1, and c- 
FOS.l. Differences in patient demographics, pathology 
and treatment, non-standard tissue handling, experimen- 
tal and statistical methods may all contribute to differ- 
ences in gene lists. Differences with other published gene 
lists might also reflect the fart that in our study, only sam- 
ples from patients with androgen-insensitive prostate can- 
cer were used. Additionally, in our experimental design 
we have incorporated features that increase the signifi- 
cance of our findings and increase the likelihood that the 
genes identified truly reflect the biology of metastatic 
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fable 3: Transcripts with median values with at least 3 fold difference between metastatic and primary tumor samples 



.Gene Symbol 


ProbeJD 


PI 


P2 


P3 


P4 


Upregulated Transcripts 














22.37 


5.78 


13.25 


56.28 


sm 




24.16 


26.78 


4.75 


5.39 


HBAI///HBA2 


3l525~7at 


15.14 


4.47 


13.65 


108.11 


1GR4 


43585_at 


7.39 


7.43 


20.89 


24.82 






14.35 


12.97 


12.24 


14.78 


PRO 1 073 


49666 s at 


4.56 


13.25 


10.01 


13.5 






10.1! 


6.02 


12.11 


16.31 




59070"" 


5.37 


9.69 


12.08 


13.67 


;SDCCAG3 


430l4~at 


7.99 


8.57 


11.24 


17.32 


LOC64744 


42739_at 


7.3 


9.64 


9.57 


14.51 






5.06 


4.14 


12.12 


22.01 


Sppi 


2092~s~at 


14.05 


12.94 


3.35 


4.07 


UBE2H 


5877 ~ S ~^ 


9.5 


7.45 


6.55 


15.13 


SRPKi 


63687 at 


6.06 


4.36 


10.61 


12.82 


NCK2 




5 


9.14 


8.5 


7.34 


HISTI H3H 




7.26 


17.07 


8.47 


5.61 


PPP4R2 


Jt 


5.09 


6.97 


8.59 


16.16 


p8orfl6 


~ at 


6.54 


8.15 


9.53 


7.39 




~ at 


3.41 


7.47 


15.31 


8.01 






8.25 


6.42 


7.04 


10.6 


Ip400 




5.94 


9.27 


4.51 


9.32 


GOLTIA 


45 1 44_at 


3.9 


6.17 


8.37 


12.32 




44791^" 


9.83 


7.1 


6.25 


7.03 


LOC2840S8 




8.25 


10.17 


4.55 


5.86 




580 ^ 


3.42 


6.04 


8.03 


11.32 


NFATC2IP 


38864at 


3.26 


4.83 


9.58 


9.19 


SELIL 


40689 at 


4.71 


7.84 


6.13 


10.94 


HLLT2 9 




3.43 


6.26 


8.92 


7.52 




65205~at 


3.43 


7.13 


6.57 


13.01 


SC4MOL 


46802~at 


22.91 


7.35 


5.62 


6.17 




62671 at 


6.38 


7.13 


5.74 


11.03 


&R.C6 


46558 _ at 
_a 


5.67 


8.59 


7.5 


5.92 


MAP4K4 


~ at 


4.86 


4.32 


8.7 


8.52 


MLLT2 


~ 3t 


4.65 


3.99 


9.39 


8.1 


k 


Jt 


8.71 


5.94 


6.35 


6.25 




4.23 


4.87 


7.39 


8.23 


ACAS2 


62783~at 


4.29 


6.14 


7.02 


5.9 


1. 


60658_at 


3.4 


6.84 


5.19 


9.1 


SUMOI 


4955 l_at 


4.05 


7.2 


4.77 


7.75 


AR 


I578_g_at 


7.56 


4.86 


5.32 


6.37 


GALNT7 


59IOI_at 


8.41 


4.12 


5.11 


6.54 


(3PR75 


44203_at 


5.14 


8.32 


3.9 


6.31 


TBLIXRI 


6500l_r_at 


3.53 


12.3 


4.06 


7.19 


HSDI7BI2 


43292_at 


4.74 


8.88 


3.64 


6.28 


l*IRPS28 


43095_at 


5.79 


5.39 


5.58 


5.14 


FN 1 


647l9_at 


27.07 


6.02 


4.05 


4.93 


GPRI58 


442l4_at 


7.21 


3.33 


4.32 


6.62 




48069_at 


6.27 


9.88 


3.38 


4.55 


FLJ2I657 


58778_at 


4.34 


5.6 


6.17 


5.18 


MLL5 


4330 l_at 


4.76 


3.61 


5.87 


10.34 




5576 l_at 


3.78 


4.88 


5.65 


6.93 
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table 3: Transcripts with median values with at least 3 fold difference between metastatic and primary tumor samples (Continued) 



;DLGI 


4723 1 _at 


3.4 


4.77 


6.22 


5.7 


MY05B 


6328l_r_at 


3.29 


6.I7 


4.29 


6.84 




49268_at 


3.55 


1 9.86 


3.61 


6.75 


FUS 


4350l_at 


3.93 


3.78 


6.42 


8.97 


CCDC35 


54684_at 


4.9 


8.I4 


3.55 


5.43 




4343S_at 


6.85 


4.83 


4.82 


5.49 


;SMA4 


3292 1 _at 


4.68 


5.53 


5.74 


4.26 


•NCOAI 


459S3_at 


6.53 


4.I3 


3.58 


6.06 


;S I O0A8 


4l096_at 


4.22 


5.89 


3.8 


22.58 


IPRKCBPI 


53493_at 


4.65 


7.37 


4-5 


5.35 




6S083_at 


3.I8 


3.96 


6.01 


9.19 


CAMSAPI 


62630_at 


4.45 


5.8 


3.36 


5.36 


EEFIG 


4l903_at 


5.I9 


4.58 


4.31 


5.34 


,'EIFS 


5l379_at 


3.44 


4.08 


5.62 


11.07 


MAML3 


49879_at 


3.39 


3.22 


10.27 


5.87 


C2lorfl06 


5965 1 _at 


3.I9 


4.02 


5.23 


6.44 


VCIPI35 


427l5_at 


3.37 


3.6 1 


5.52 


8.55 


FOX03A 


55502_at 


3.48 


4.37 


6.97 


4.74 




49l43_s_at 


4.23 


4.62 


4.41 


5.78 


CNMT 


46482_at 


3.59 


4.84 


4.24 


4.64 


DONSON 


48549_at 


4. 1 


3.58 


4.66 


5.28 




43436_g_at 


4.98 


3.75 


3.58 


5.09 


PKP4 


66327_at 


3.3 1 


3.88 


4.56 


6.2 


PCBP2 


55393_at 


3.73 


3.I9 


4.36 


6.29 


jbPEB4 


57l69_at 


3.7 


3.92 


4.14 


4.48 


CUGBPI 


34683_at 


4.26 


3.76 


3.13 


4.78 


IrALZ 


47458_at 


4.2 1 


3.65 


3.82 


4.09 




5l586_at 


3.5 1 


4 


4.99 


3.89 


RALA 


39253_s_at 


3.92 


4.3 


3.29 


3.85 


MLL5 


45092_at 


4.36 


3.2I 


4.48 


3.39 


PABPCI 


44806_at 


3.74 


3.98 






EIFIAX 


34278_at 


3.99 


3.47 


3.84 


3.19 


C7orf2 


42l73_at 


3.I5 


3.27 


5.07 


4.04 




63l47_at 


3.25 


5.4 


3.12 


4.04 


RAD23B 


4 1 1 57_at 


3.2 


3.46 


3.64 


4.45 




6l037_at 


3.44 


3.56 


3.47 


3.73 


NFATC I 


39l43_at 


3.I3 


3.2 1 


9.06 


3.78 


JARIDIA 


50532_at 


3.22 


3.32 


3.54 


4.12 


PDUM5 


37366_at 


3.02 


3.58 


3.42 


3.16 


Downregulated Ti 


ranscripts 










NEFH 


33767_at 


-I I7.I5 


-1 47.36 


-9.9 


-17.18 


ClOorfl 16 


32527_at 


-35.49 


-29.63 


-46.85 


-66.5 


KLK \ I 


40035_at 


-23.65 


-1 9.24 


-39.73 


-62.15 


FAM3B 


59657_at 


-I5.8I 


-27.92 


-26.09 


-25.97 


PGM5 


52l40_at 


-23.87 


-26.5 


-44.27 


-17.72 


MRGPRF 


52946_at 


-I5.6I 


-1 8.57 


-30.59 


-70.95 


l<RTI5 


37582_at 


-2 1. 85 


-20.74 


-19.22 


-33.68 


PTN 


34820_at 


-II.62 


-3 1 .95 


-10.24 


-27.11 


SELM 


64449_at 


-6.36 


-8.4 


-29.23 


-39.36 


Mylk 


46276_at 


-5.87 


-1 5.22 


-22.57 


-20.86 


5YNP02 


5036l_at 


-I5.I4 


-1 5.77 


-20.15 


-84.14 


(CRT5 


6l3_at 


-I3.2I 


-II. 12 


-22.66 


-32.96 


FOS 


2094_s_at 


-1 0.72 


-25.75 


-13.72 


-16.45 


PKPI 


5l2l4_at 


-II.57 


-16.34 


-11.83 


-17.85 



Page 11 of 21 
(page number not for citation purposes) 



BMC Cancer2QQ7, 7:64 http://www.biomedcentral.eom/1471-2407/7/64 



Table 3: Transcript 




with at least 3 fold differ 


ence between metastatic 


: and primary tumoi 


- samples (Continued) 




4292 l_at 


-9.96 


-1 1.67 


-15.61 


-16.5 


|RAB34 


45269_at 


-14.36 


-1 1.54 


-17.49 


-10.35 


ir~ 


48927_at 


-10.61 


-14.93 


-8.77 


-21.91 


ALOXI5B 


37430_at 


-12.47 


-12.41 


-14.17 


-9.1 


FOS 


l9IS_s_at 


-7.59 


-26.38 


-1 1.03 


-12.1 1 


TMEMI6G 


62387_at 


-9.63 


-13.32 


-12.59 


-9.93 


J" 


64676_at 


-17.3 


-9.39 


-6.32 


-13.05 


5FRPI 


3252 l_at 


-13.1 


-5.73 


-8.29 


-16.73 


NDFIP2 


605IO_at 


-7.2 


-9.23 


-11.72 


-15.15 


FHOD3 


50298_at 


-9.96 


-12.84 


-5.59 


-10.96 


WNT5B 


6i292_s_at 


-8.72 


-1 1.85 


-5.42 


-13.92 


SYNP02 


48039_at 


-1 1.04 


-8.8 


-12.64 


-9.34 


BOC 


64423_s_at 


-3.63 


-8.16 


-1 1.8 


-54.66 


&LC20A2 


I I37_at 


-9.27 


-5.08 


-10.51 


-12.61 


COL8A2 


52652_g_at 


-7.95 


-9.99 


-11.56 


-9.75 




52678_at 


-9.69 


-9.99 


-3.76 


-17.93 


FOS 


I9l6_s_at 


-7.58 


-21.81 


-6.93 


-1 1.58 


ARGBP2 


SI939_at 


-7.77 


-13.86 


-10.4 


-8.71 


CTGF 


64342_at 


-4.21 


-4.15 


-20.44 


-14.87 


EPHB6 


39930_at 


-8.61 


-9.66 


-8.32 


-19.41 


SYNP02 


60532_at 


-9.77 


-5.54 


-8.77 


-9.03 


iNR4AI 


280_g_at 


-8.68 


-13.49 


-5.82 


-8.58 


PKFZPS64O0823 


54033_at 


-4.67 


-3.72 


-11.83 


-20 


GST02 

F" 


4S609_at 


-4.73 


-6.81 


-9.6 


-16.18 


4932 l_at 


-7.91 


-8.41 


-9.24 


-3.88 


EGR3 


40375_at 


-9.89 


-7.71 


-8.49 


-6.44 


SYNP02 


6l68l_at 


-7.85 


-8.33 


-4.56 


-18.57 


PII5 


5836 l_at 


-3.59 


-4.26 


-12.77 


-11.74 


fOSB 


36669_at 


-8.81 


-6.27 


-7.6 


-8.39 


QGN 


43507_g_at 


-3.56 


-8.26 


-7.19 


-25.54 


MOXDI 


36834_at 


-5.4 


-1 1.7 


-10 


-3.85 


LSAMP 


43930_at 


-3.05 


-7.62 


-9.76 


-7.67 


EGR2 


37863_at 


-7.7 


-5.52 


-7.23 


-15.41 


PKFZp686D0853 


49770_at 


-10.18 


-7.66 


-7.16 


-4.39 


L.GPI 


52826_at 


-13.75 


-5.94 


-3.83 


-8.1 1 


ME3 


352l6_at 


-7.45 


-9.26 


-6.54 


-5.32 


PPPIRI4A 


58774_at 


-6.68 


-6.14 


-7.31 


-7.87 


FLJ22386 


50l98_at 


-6.8 


-3.64 


-6.98 


-6.65 


NR4AI 


279_at 


-5.31 


-8.04 


-5.1 1 


-8.48 


WFDCI 


64lll_at 


-3.79 


-1 1.21 


-6.64 


-6.66 


ZFP36 


40448_at 


-6.39 


-6.86 


-7.25 


-3.61 


GACHDI 


43554_at 


-6.68 


-3.34 


-17.46 


-6.57 


RLNI 


35070_at 


-6.78 


-1 1.78 


-5.14 


-6.39 


f~ 


49975_at 


-6.43 


-6.16 


-6.74 


-10.1 1 


tYBRDI 


65852_at 


-6.43 


-4.79 


-6.7 


-7.23 


PER3 


53766_at 


-15.43 


-6.79 


-5.56 


-6.29 


|lNI 


37283_at 


-4.47 


-7.36 


-5.55 


-7.48 


DNCI2 


35788_at 


-4.2 


-8.68 


-3.02 


-10.64 


MRVII 


43966_at 


-6.76 


-5.28 


-12.19 


-6.09 


AZGPI 


35834_at 


-6.32 


-3.86 


-38.18 


-6.18 


MGC 14839 


48949_at 


-8.96 


-4.19 


-8.25 


-3.61 


SMTN 


64499_s_at 


-5.2 


-15.22 


-7.18 


-4.42 


HSPCI57 


50l79_at 


-5.66 


-3.18 


-6.63 


-8.09 


WFDC2 


33933_at 


-5.3 


-6.5 


-5.73 


-6.81 


BJG2 


36634_at 


-6.99 


-3.13 


-9.25 


-5.22 
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fable 3: Transcripts with median values with at least 3 fold difference between metastatic and primary tumor samples (Continued) 



AXIN2 


64l29_at 




-4.97 


-6.97 


-7.18 


-4.2 


PDGFC 


452l7_at 




-4.32 


-7.53 


-8.81 


-3.97 


MULTIO 


63345_at 




-7.2 


-5.85 


-5.9 


-3.84 


BMP7 


49273_g_at 




-4.58 


-4.89 


-6.82 


-13.13 


MCC 


49504_r_at 




-5.9 


-5.71 


-5.08 


-5.84 








-8.15 


-5.65 


-4.18 


-5.88 


GSTT2 


mwVat 




-6.47 


-5.05 


-6.8 


-4.66 


sspN , " 


6S647_at 




-5.4 


-5.88 


-3.12 


-17.61 


upK3A 






-5.37 


-4.71 


-5.81 


-6.91 


L DESA 


54668"" 




-4.44 


-5.17 


-5.87 


-9.56 


>SD3 


63832~at 




-3.19 


-6.04 


-4.98 


-6.58 


ALDH7AI 


6| „ ~ a 




-5.85 


-5.14 


-5.88 


-3.13 


FMOD 


3343 l~at 






-4.3 


-4.9 


-6.04 


TSPAN2 


53693~at 




-6.38 


-4.49 


-4.54 


-6.77 


DKFZP586H2I23 


40017 at 




-6.52 


-6.49 


-4.32 


-3.91 




33883_at 




-5.43 


-3.58 


-6.35 


-5.18 


PODN 






-4.16 


-5.3 


-4.84 


-4.98 


buspi 






-6.53 


- 1 6.66 


-3.02 


-3.19 


SLC22AI7 


5889 ~ a s a( . 




-4.93 


-5.81 


-4.66 


-4.44 


CDHIO 


47535~at 




-4.87 


-3.19 


-8.27 


-4.65 




64I63 _ at 




-3.66 


-5.03 


-4.79 
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Table 4: GO and pathway annotation of genes and pathways altered in rr 
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Gene Symbol 
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*Bone remodeling 














j34342_s_at 


SPPI 


24.7 


24.16 
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-2.5 


-5.33 
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3.36 


60S68_at 


BCL2 


2.07 


2.18 
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2.5 
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2.53 
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-25.75 


-13.72 


-16.45 



Genes whose median value shows at least 2 fold change in every metastatic patient were annotated by NIH's DAVID tool. Genes belonging to 
Selected GO categories with a significant number of differentially expressed genes are shown. Overa//_FC, the fold change in the mean expression 
value of all metastatic samples compared to the primary samples;FC_P/, FC_P2, FC_P3, FC_P4; a patient-specific median expression value was 
calculated using all the samples from each patient This was compared to the median value for the primary tumors. Patient P4 had only sample and 
this was considered the median. 
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Figure 3 

Comparison of Z-transformed expression values between the Affymetrix and Codelink platforms. Gene expres- 
sion data from Affymetrix and CodeLink experiments was Z-transformed to allow comparison. Data for selected differentially 
expression genes is shown. 
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Figure 4 

Validation of differentially expressed genes with quantitative real time PCR. QPCR was performed on RNA from 
samples used in the CodeLink analysis. Three selected metastatic RNA samples from each patient were pooled together 
(except for patient FB666 n = I) and therefore four RNA samples, each representing one metastatic patient and 5 primary 
tumors were tested. The insert shows average expression values for metastatic and primary tumors. 



prostate cancer. First, by subtracting transcripts previously 
identified as being expressed by the prostatic stroma, we 
have incorporated previous knowledge about the expres- 
sion profiles of different components of prostate tumors 
in order to focus on those transcripts intrinsic to meta- 
static cells. This takes into account the fact that metastatic 
tumors do not contain all the tissue elements present in 
o'rgan-confined tumors. The major benefit of this strategy 
is to better define the genes that are down-regulated in 
metastatic tumor cells. Second, by analyzing multiple 
tumor samples from each patient, we have addressed the 
fact that metastatic prostate cancer shows significant het- 
erogeneity, even within the same patient [21]. This is cor- 
roborated by our results, and we address this issue by 



focusing our analysis on the transcripts that show signifi- 
cant differential expression in all metastatic sites within 
and between patients. 

Multiple biological processes appear to be altered in met- 
astatic prostate cancer. One common theme that has 
emerged from studies of metastatic disease is the central 
role of the androgen receptor in the development of 
androgen resistant disease. Several mechanisms including 
amplification of the AR gene, upregulation of mRNA 
expression to allow binding by low levels of androgens, 
mutations in the ligand binding domain (LBD) that allow 
the receptor to be activated by antagonists, and alteration 
in the normal AR signaling pathway, have been proposed 
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to explain the ability of prostate cancer to recur in the 
presence of androgen ablation therapy [33,34]. Consist- 
ent with previous observations, AR is up-regulated in all 
metastatic samples in our study. Similarly, gene expres- 
sion changes of the MAP kinase pathway in metastasis 
may be related to the development of the AARPC pheno- 
type. The gene list from our analysis shares some similari- 
ties with mouse xenograft prostate cancers models 
(CWR22) of androgen independence. MSMB, CCND1, 
EFNA3, FKBP and ADM, HGF ate similarly regulated in 
mouse models and in our study [35-37]. 

Changes in the expression level of several additional tran- 
scripts may reveal clues about the mechanism of metasta- 
sis and androgen resistance. We find upregulation of the 
enzyme UGT2B15 in all metastatic patients. Upregulation 
of UGT2B15 in androgen independent prostate cancer has 
been reported previously [14]. This increase appears para- 
doxical, since UGT2B15 is involved in hormone inactiva- 
tion. However, as suggested by Stanbrbugh et dl [14], 
upregulation of multiple genes related with androgen 
metabolism might reflect that metastatic tumor cells have 
an increased capacity to convert weak androgens into tes- 
tosterone or DHT. However, in contrast to their findings, 
transcripts for AKR1C3, SRD5A1, HSD3B2, AKR1C2, 
AtfCRl CI are not consistently upregulated in the metastatic 
samples in our study. Interestingly, when reviewing indi- 
vidual values for each sample, some metastases indeed 
show higher levels for some of these transcripts, which 
niight reflect the heterogeneity of metastatic prostate can- 
cer phenotypes. Another possible explanation for this dis- 
crepancy is that the metastatic samples used by 
Stanbrough et dl. are all from bone metastases and this 
type of sample is not represented in our study. Clearly, 
further investigation into the role of these pathway genes 
in the development of androgen resistance by metastatic 
samples is needed. 

Several genes involved in cell-cell interaction and cell 
aidhesion appear to be up-regulated in these tumors. SPP1 
(psteopontin), a secreted, integrin-binding glycoprotein 
With adhesive properties, has been shown to be correlated 
with metastasis to the bone and with poor prognosis in 
various cancers and is highly upregulated in all the meta- 
static samples in our study. Elevated plasma osteopontin 
levels have also been correlated with lower survival and 
bone metastasis in hormone resistant prostate cancer 
[25]. Interestingly, Stanborough and collaborators also 
identified SPP1 as upregulated in their metastatic sam- 
ples, however, their interpretation was that this increase 
w;as part of the bone response to the metastases. Our study 
confirms that upregulation of SPP1 is a feature intrinsic to 
androgen-resistant metastatic prostate cancer, independ- 
ent of the site of metastasis. It has been postulated that 
metastasis to specific target organs may require not only 



expression of SPP1 but an additional set of signaling mol- 
ecules that promote metastasis to the specific organ. SPP1 
when expressed with IL11 has been shown to promote 
metastasis of breast cancer cells to the bone [38] but not 
to the adrenal medulla. Further detailed studies are 
required to address the specific role of SPP1 and other co- 
expressed genes in prostate cancer metastasis and whether 
SPP1 represents a potential therapeutic target for andro- 
gen-resistant disease. Interestingly, the gene expression 
profile termed as "bone module" and postulated as a hall- 
mark of tumor metastasis to bone [39] is not dysregulated 
in our study, most likely reflecting the fart that we did not 
assay bone metastatic samples. It is also possible, that the 
role that SPP1 plays in metastasis to bone and/or other 
organs may involve distinct mechanisms [38]. 

Metastatic tumors have been described as undergoing an 
epithelial to mesenchymal transition with loss of the dif- 
ferentiated phenotype. Downregulation of transcription 
factors such as JUN has been observed in advanced stages 
of other cancers and its loss of activity has been postulated 
to be involved in this transition [40]. In our study, both 
FOS and JUNB, which are upregulated in primary tumors 
compared to normal prostate tissue are highly downregu- 
lated in the metastatic samples. FBN, also representative 
of the EMT transition [38] is overexpressed in our meta- 
static samples. Our analysis has also identified a number 
of additional genes, such as KLK1 1,STC1 and S100A8 that 
are uniformly regulated in all metastatic patients. The role 
of S100A8 in prostate cancer has been studied with evi- 
dence suggesting that it is elevated in prostate cancer and 
may be involved in MAP kinase and NFK-B signalling 
[41,42], STC1, involved in calcium homeostasis, has been 
reported to have osteoblastic and angiogenic modulator 
properties with altered expression in some cancers [43- 
45]. The serine protease KLK1 1 appears to be regulated in 
prostate cancer with negative correlation between aggres- 
siveness and expression [46]. 

A recent study observed overexpression of 62 genes due to 
surgical manipulation related ischemia of the prostate 
[47], In our study, 12 out of the 62-gene ischemia profile 
are downregulated in all metastatic samples. This gene list 
includes DUSP1, BTG2, IER2, PTGS2, NR4A1, AMD1, 
C20orf35, KLF4, RAB4A, KLF6, CTGF and GOLPH2. In our 
data set, these genes represent only 0.01% of the total 
number of genes differentially regulated in all metastatic 
samples. Since our metastatic samples all originate from 
autopsy studies, it is likely that they had been exposed to 
longer ischemia than the organ confined samples 
obtained from surgical specimens. Thus, if the differences 
we observed were related to the ischemia, we would have 
expected an increase in the expression of these genes, and 
not the observed downregulation. Therefore, it is unlikely 
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that surgical manipulation can explain the differential 
gene expression between metastatic and primary tumors. 

Conclusion 

In. summary, our results support the roles for specific cell 
adhesion, androgen metabolism and transcription factor 
genes in the development of androgen-independent met- 
astatic prostate cancer. Furthermore, the differentially 
expressed transcripts in metastatic tumors that we report 
have been validated with two independent sets of primary 
tumors, two gene expression microarray platforms, and 
selected genes were further validated by qRT-PCR. Our 
results corroborate the notion that metastatic prostate 
cancer is quite heterogeneous within a single patient. 
Despite this heterogeneity our experimental design 
allowed us to identify common expression profiles for 
androgen-independent metastatic prostate cancer. 
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